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Abstract: We consider AA = 1 supersymmetric U(N) field theories in four dimensions with adjoint 
chiral matter and a multi-trace tree-level superpotential. We show that the computation of the 
effective action as a function of the glueball superfield localizes to computing matrix integrals. 
Unlike the single-trace case, holomorphy and symmetries do not forbid non-planar contributions. 
Nevertheless, only a special subset of the planar diagrams contributes to the exact result. Some of 
the data of this subset can be computed from the large- A limit of an associated multi-trace Matrix 
model. However, the prescription differs in important respects from that of Dijkgraaf and Vafa 
for single-trace superpotentials in that the field theory effective action is not the derivative of a 
multi-trace matrix model free energy. The basic subtlety involves the correct identification of the 
field theory glueball as a variable in the Matrix model, as we show via an auxiliary construction 
involving a single-trace matrix model with additional singlet fields which are integrated out to 
compute the multi-trace results. Along the way we also describe a general technique for computing 
the large- N limits of multi-trace Matrix models and raise the challenge of finding the field theories 
whose effective actions they may compute. Since our models can be treated as M = 1 deformations 
of pure N = 2 gauge theory, we show that the effective superpotential that we compute also follows 
from the N = 2 Seiberg-Witten solution. Finally, we observe an interesting connection between 
multi-trace local theories and non-local field theory. 
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1. Introduction 

Dijkgraaf and Vafa have recently made the remarkable proposal that the superpotential and other 
holomorphic data of M = 1 supersymmetric gauge theories in four dimensions can be computed 
from an auxiliary Matrix model ||, ^ ]3|. While the original proposal arose from consideration 
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of stringy dualities arising in context of geometrically engineered field theories, two recent papers 
have suggested direct field theory proofs of the proposal j|, |5|. These works considered U(N) 
gauge theories with an adjoint chiral matter multiplet $ and a tree-level superpotential W(3?) = 
g& Tr(<l? fc ). Using somewhat different techniques (Q uses properties of superspace perturbation 
theory while || relies on factorization of chiral correlation functions, symmetries, and the Konishi 
anomaly) these papers conclude that: 

1. The computation of the effective superpotential as a function of the glueball superfield reduces 
to computing matrix integrals. 

2. Because of holomorphy and symmetries (or properties of superspace perturbation theory), 
only planar Feynman diagrams contribute. 

3. These diagrams can be summed up by the large-iV limit of an auxiliary Matrix model. The 
field theory effective action is obtained as a derivative of the Matrix model free energy. 

Various generalizations and extensions of these ideas (e.g., M = 1* theories |], [7j, fundamental 
matter Js], |J, quantum moduli spaces Juj, non-super symmetric cases [11], other gauge groups 
[12, [b^, [u], baryonic matter If, 17], gravitational corrections [18, 19], and Seiberg Duality 



1 20, ^]) have been considered in the recent literature. 

A stringent and simple test of the Dijkgraaf-Vafa proposal and of the proofs presented in [||, || 
is to consider superpotentials containing multi-trace terms such as 

W($) = 52 Tr($ 2 ) + 54 Tr($ 4 ) + g 2 (Tr($ 2 )) 2 . (1.1) 

We will show that for such multi-trace theories: 

1. The computation of the effective superpotential as a function of the glueball superfield still 
reduces to computing matrix integrals. 

2. Holomorphy and symmetries do not forbid non-planar contributions; nevertheless only a 
certain subset of the planar diagrams contributes to the effective superpotential. 

3. This subclass of planar graphs also contributes to the large- N limit of an associated multi- 
trace Matrix model. However, because of differences in combinatorial factors, the field theory 
effective superpotential cannot be obtained simply as a derivative of the multi-trace Matrix 
model free energy as in ||. 

4. Multi-trace theories can be linearized in traces by the addition of auxiliary singlet fields A4. 
The superpotentials for these theories as a function of both the A{ and the glueball can be 
computed from an associated Matrix model. This shows that the basic subtlety involves the 
correct identification of the field theory glueball as a variable in a related Matrix model. 

The plan of this paper is as follows. In Sec. 2 we carefully analyze the methods of Q and 
generalize them so that they apply to an J\f = 1 U (N) gauge theory in four dimensions with a tree- 



level superpotential of the form (1.1). Along the way we introduce some new techniques that deepen 
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our understanding of the selection rules determining which perturbative field theory diagrams 
contribute to the effective superpotential of an N = 1 field theory. Using this understanding we 
demonstrate how the conclusions (1) and (2) above arise and show that contributing diagrams are 
tree-like graphs in which single-trace diagrams are pasted together by double-trace vertices through 
which no momentum flows. We illustrate our results by explicitly computing the superpotential to 
the first few orders in perturbation theory. Finally, we observe an intriguing connection between 
multi-trace local theories and non-local field theory. 

Since the field content of pure N = 2 supersymmetric U (N) gauge theory in four dimensions 
consists in M = 1 language of a vector multiplet W a and an adjoint chiral multiplet the 
superpotential can be treated as a deformation of an j\[ = 2 theory to an j\[ = 1 theory. 

Hence, we can use global symmetries, holomorphy, regularity conditions, and the Seiberg-Witten 
solution of M = 2 gauge theory to compute the exact superpotential. We carry out this procedure in 
Sec. 3, using the fact that the vacuum expectation value of the product of chiral operators (Tr($ 2 ) 2 ) 
factorizes as (Tr(<3? 2 )) 2 . We show that the result exactly captures the subset of the planar diagrams 
that contribute to the exact field theory superpotential. The assumption of factorization in the 
Seiberg-Witten analysis is equivalent to the vanishing of a certain subset of planar diagrams in our 
perturbative computations. 

In Sec. 4 we demonstrate a general technique for solving U(M) matrix models (or general 
complex matrix models) with multi-trace potentials. The essential observation, following Das, 
Dhar, Sengupta, and Wadia 22], is that in the large-M limit, mean field methods can be used 
to solve for the effect on a single matrix eigenvalue of the rest of the matrix. We explain the 
general method and solve two examples in detail. The first example has a potential V(<&) = 
M(5 2 Tr($ 2 ) + 5f 4 Tr($ 4 ) + |fTr($ 2 ) 2 ) for (j) e U(M). By expanding the exact large-M result in 
powers of the couplings we demonstrate how this limit computes the data relevant for a certain 
subset of the planar contributions to the effective action of the field theory with the tree-level 
superpotential in ( |1 . 1|) . In the proposal of Dijkgraaf and Vafa || and the subsequent generalizations 
(e.g., H to the field theory effective action was related simply to the free energies of auxiliary 
matrix models and their derivatives. We demonstrate the absence of such a relation between multi- 
trace field theories and multi-trace Matrix models. As a further illustration of the mean field 
technique for computing large-M limits, we study a matrix model with a general quartic potential. 

Multi-trace field theories can be linearized in traces by the introduction of new singlet fields 
which can be integrated out to produce the multi-trace theory. In Sec. 5 we show how this procedure 
is carried out and relate the resulting linearized superpotential to a Matrix model following the 
techniques of 0]. In the Matrix model integrating out the singlets at the level of the free energy 
reproduces the multi-trace results that do not agree with the field theory. However, integrating out 
after computing the linearized field theory superpotential leads to agreement. This shows that the 
basic subtlety here involves the correct identification of the field theory glueball as a variable in an 
associated Matrix moedl. 

It is worth mentioning several further reasons why multi-trace superpotentials are interesting. 
First of all, the general deformation of a pure N = 2 field theory to an J\f = 1 theory with ad- 
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joint matter involves multi-trace superpotentials, and therefore these deformations are important 
to understand. What is more, multi-trace superpotentials cannot be geometrically engineered |23|] 
in the usual manner for a simple reason: in geometric engineering of gauge theories the tree-level 
superpotential arises from a disc diagram for open strings on a D-brane and these, having only one 
boundary, produce single-trace terms. In this context, even if multi-trace terms could be produced 
by quantum corrections, their coefficients would be determined by the tree-level couplings and 
would not be freely tunable. Hence comparison of the low-energy physics arising from multi-trace 
superpotentials with the corresponding Matrix model calculations is a useful probe of the extent 
to which the Dijkgraaf-Vafa proposal is tied to its geometric and D-brane origins. In addition to 
these motivations, it is worth recalling that the double scaling limit of the U(N) matrix model 
with a double-trace potential is related to a theory of two-dimensional gravity with a cosmological 
constant. This matrix model also displays phase transitions between smooth, branched polymer 
and intermediate phases |p!f |. It would be interesting to understand whether and how these phe- 
nomena manifest themselves as effects in a four dimensional field theory. The results of our paper 
suggest that these phase transitions and the physics of two-dimensional cosmological constants are 
embedded within four-dimensional field theory. It would be interesting to explore this. Finally, 
multi-trace deformations of field theories have recently made an appearance in the contexts of the 
AdS /CFT correspondence and a proposed definition of string theories with a nonlocal worldsheet 
theory pi. 



2. Multi- Trace Superpotentials from Perturbation Theory 

In this section we begin by reviewing the field theoretic proof that when treated as a function 
of the glueball superfield, the effective superpotential of an N = 1 supersymmetric gauge theory 
with single-trace tree- level interactions is computed by planar matrix diagrams ||. We will 
then describe how these arguments are modified by the presence of multi-trace terms in the tree- 
level action. Finally, we will explicitly illustrate our reasoning by perturbatively computing the 
diagrams that contribute to the effective superpotential of a multi-trace theory up to third order 
in the couplings. We will always work around a vacuum with unbroken U(N) symmetry. 

2.1 A Schematic Review of the Field Theory Superpotential Computation 

Below we give a schematic description of the methods of 0] for the computation of the effective 
superpotential of an J\f = 1 field theory. While discussed theories with single-trace Lagrangians, 
we will find that most of their arguments will generalize easily to multi-trace theories. 

1. The Action: The matter action for an J\f = 1 U{N) gauge theory with a vector multiplet V, 
a massive chiral superfield and superpotential W(<&), is given in superspace by 

£($,$) = [ d 4 xd 4 e$e v $+ fd 4 xd 2 6W($) + h.c. (2.1) 
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2. The Goal: We seek to compute the effective superpotential as a function of the glueball 

superfield 

S = ^Tr(W a W a ), (2.2) 

where 

W a = W 2 e- v D a e v (2.3) 

is the gauge field strength of V, with D a = d/d6 a and = d/d6 a +i6 a d a a the superspace 
covariant derivatives, and D 2 = ^D a D a and D 2 = ^D a Da- The gluino condensate S is a 
commuting field constructed out of a pair of fermionic operators W a . 

3. The Power of Holomorphy: We are interested in expressing the effective superpotential in 

terms of the chiral glueball superfield S. Holomorphy tells us that it will be independent 
of the parameters of the anti-holomorphic part of the tree-level superpotential. Therefore, 
without loss of generality, we can choose a particularly simple form for W($): 

= ^m<l 2 . (2.4) 

Integrating out the anti-holomorphic fields and performing standard superspace manipula- 
tions as discussed in Sec. 2 of Q, gives 

S = f d 4 xd 2 6 f-y=* (□ - iW a D a ) $ + WtreeW) (2-5) 

as the part of the action that is relevant for computing the effective potential as a function 
of S. Here, □ = \d a o t d aa is the d'Alembertian, and Wt re e is the tree-level superpotential, 
expanded as ^m$ 2 + interactions. (The reader may consult Sec. 2 of for a discussion of 
various subtleties such as why the □ can be taken as the ordinary d'Alembertian as opposed 
to a gauge covariantized Dcov)- 

4. The Propagator: After reduction into the form fl2.5|) , the quadratic part gives the propaga- 

tor. We write the covariant derivative in terms of Grassmann momentum variables 



D a = d/d6 a := -in a , (2.6) 

and it has been shown in 0] that by rescaling the momenta we can put m = 1 since all 
m dependence cancels out. Then the momentum space representation of the propagator is 
simply 

poo 

/ dsi exp (-Si(p 2 + W a ir ia + m)) , (2.7) 
J o 

where is the Schwinger time parameter of i-th Feynman propagator. Here the precise form 
of the W a ir a depends on the representation of the gauge group that is carried by the field 
propagating in the loop. 
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5. Calculation of Feynman Diagrams: The effective superpotential as a function of the glue- 
ball S is a sum of vacuum Feynman diagrams computed in the background of a fixed constant 
W a leading to insertions of this field along propagators. In general there will be I momentum 
loops, and the corresponding momenta must be integrated over yielding the contribution 



J] d 4 p a dsi e~ s ^ I • ( / II d ^ ds * e ~ 




— Iboson ' Ifermion ' p C^"^) 
777/ 

to the overall amplitude. Here a labels momentum loops, while i = 1, . . . , P labels propaga- 
tors. The momenta in the propagators are linear combinations of the loop momenta because 
of momentum conservation. 

6. Bosonic Momentum Integrations: The bosonic contribution can be expressed as 

74 



/dp 
n ^ exp 



^2PaM ab (s)p b 



a.b 



(2.9, 



(4vr)^ (det M{s)f 



where we have defined the momentum of the 7-th propagator in terms of the independent 
loop momenta p a 

Vi = Y^ L ^ ( 2 - 10 ) 

a 

via the matrix elements L{ a £ {0, ±1} and 

M ab {s) =^2 Si L ia L ib . (2.11) 

i 

7. Which Diagrams Contribute: Since each momentum loop comes with two fermionic n a 

integrations ( |2.8[ ) a non-zero amplitude will require the insertion of 2£ 7r a s. From fl2.7| ) we 
see that that ttq, insertions arise from the power series expansion of the fermionic part of the 
propagator and that each ir a is accompanied by a W a . So in total we expect an amplitude 
containing 21 factors of W a . Furthermore, since we wish to compute the superpotential as a 
function of S ~ Tr(W a W a ) each index loop can only have zero or two W a insertions. These 
considerations together imply that if a diagram contributes to the effective superpotential as 
function of the S, then number of index loops h must be greater than or equal to the number 
of momentum loops £, i.e., 

h>£. (2.12) 

8. Planarity: The above considerations are completely general. Now let us specialize to U(N) 

theories with single-trace operators. A diagram with I momentum loops has 

h = t + l-2g (2.13) 

index loops, where g is the genus of the surface generated by 't Hooft double line notation. 
Combining this with (|2.12j ) tell us that g = 0, i.e., only planar diagrams contribute. 



- 6 - 



9. Doing The Fermionic Integrations: First let us discuss the combinatorial factors that 
arise from the fermionic integrations. Since the number of momentum loops is one less than 
the number of index loops, we must choose which of the latter to leave free of W a insertions. 
This gives a combinatorial factor of h, and the empty index loop gives a factor of N from the 
sum over color. For each loop with two W a insertions we get a factor of ^W a W a = 1Gtt 2 S. 
Since we are dealing with adjoint matter, the action of W a is through a commutator 

exp(- Si [Wf,-]7r ia ) (2.14) 

in the Schwinger term. (See the appendix of [^] for a nice explanation of this notation as it 
appears in Q. In Sec. 2.2 we will give an alternative discussion of the fermionic integrations 
that clarifies various points.) As in the bosonic integrals above, it is convenient to express 
the fermionic propagator momenta as sums of the independent loop momenta: 

7Fia = y j Lja^aa, (2.15) 
a 

where the Lj a are the same matrix elements as introduced above. The authors of jl also find 
it convenient to introduce auxiliary fermionic variables via the equation 

wr = J2 L ^a- (2-i6) 

a 

Here, the Lj a = ±1 denotes the left- or right-action of the commutator. In terms of the W°, 
the fermionic contribution to the amplitude can be written as 

Y,W2M ab (s)K ba 

a,b 

(2.17) 



Ifermion = Nh(16ir 2 S) e f \{d 2 K a d 2 W a exp 
= Uir) 2e NhS e (det M(s)) 2 . 



10. Localization: The Schwinger parameter dependence in the bosonic and fermionic momen- 
tum integrations cancel exactly 

Iboson ' I f ermion NHS , (2.18) 

implying that the computation of the effective superpotential as a function of the S localizes 
to summing matrix integrals. All the four-dimensional spacetime dependence has washed 
out. The full effective superpotential W e ff(S) is thus a sum over planar matrix graphs with 
the addition of the Veneziano-Yankielowicz term for the pure Yang-Mills theory [25|. The 



terms in the effective action proportional to arise exclusively from planar graphs with £ 
momentum loops giving a perturbative computation of the exact superpotential. 

11. The Matrix Model: The localization of the field theory computation to a set of planar 
matrix diagrams suggests that the sum of diagrams can be computed exactly by the large-M 
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limit of a bosonic Matrix model. (We distinguish between M, the rank of the matrices in the 
Matrix model and N, the rank of the gauge group.) The prescription of Dijkgraaf and Vafa 
does exactly this for single-trace superpotentials. Since the number of momentum loops is 
one less than the number of index loops in a planar diagram, the net result of the bosonic 
and fermionic integrations in ( |2.18| ) can be written as 

dS h 

Iboson ' Ifermion N • (2.19) 

Because of this, the perturbative part of the effective superpotential, namely the sum over 
planar diagrams in the field theory, can be written in terms of the genus zero free energy 
J~o(S) of the corresponding matrix model: 

W pert (S) = N^Fo(S), (2.20) 
T (S) =Y,^,hS h . (2.21) 

h 

This free energy is conveniently isolated by taking the large- M limit of the zero-dimensional 
one-matrix model with M x M matrices 1 3> and potential W(3>) whose partition function is 
given by 

Z - exp(M^ ) _ / m exp (-IlY *(•)) . (2.22) 

In this matrix model every index loop gives a power of M just as in the field theory compu- 
tation, and all but one index loop gives a power of S. Because of this simple fact the powers 
of the gluino condensate in the field theory superpotential can be conveniently counted by 
identifying it with the 't Hooft coupling S = Mg s , and then differentiating the matrix model 



free energy as in (2.20). Rather surprisingly the Veneziano-Yankielowicz term in W e ff(S) 



arises from the volume factor in the integration over matrices in ( p. 22 ). 



One important unanswered question is why the low-energy dynamics simplifies so much when 
written in terms of the gluino condensates. 

2.2 Computation of a Multi- Trace Superpotential 

We have reviewed above how the field theory calculation of the effective superpotential for a single- 
trace theory localizes to a matrix model computation. In this subsection we show how the argument 
is modified when the tree-level superpotential includes multi-trace terms. We consider an M = 1 
theory with the tree-level superpotential 

W tree = ^TY($ 2 ) + 24Tr($ 4 ) + £ 2 (Tr($ 2 )) 2 . (2.23) 



lr The original papers of Dijkgraaf and Vafa [[l], ^] consider M x M Hermitian matrices (i.e. matrices with real 
eigenvalues Aj). In fact, we should think of the matrices as belonging to GL(M, C) with eigenvalues distributed 
along contours in the complex plane rather than along domains on the real axis. The prior results do not depend 
crucially on this point. Indeed, they carry through exactly by analytic continuation. We thank David Berenstein for 
emphasizing this to us. 
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To set the stage for our perturbative computation of the effective superpotential we begin by 
analyzing the structure of the new diagrams introduced by the double-trace term. If 7j2 = 0, the 
connected diagrams we get are the familiar single-trace ones; we will call these primitive diagrams. 
When 7/2 7^ propagators in primitive diagrams can be spliced together by new double-trace 
vertices. It is useful to do an explicit example to see how this splicing occurs. As an example, let 




(a) (b) 

Figure 1: Two ways in which the double-trace operator: Tr(<I> 2 )Tr($ 2 ) can be contracted using the vertex 
shown in (c). 

us study the expectation value of the double-trace operator: (Tr(<I> 2 )Tr(<l> 2 )). To lowest order in 
couplings, the two ways to contract ^s give rise to the two diagrams in Figure [l]. When we draw 

these diagrams in double line notation, we find that Figure |l|a corresponding to Tr(¥$)Tr(¥$) has 

■ 1 

_ 1 1 

four index loops, while Figure |l|b corresponding to Tr(<I> <£)Tr(<I> <£) has only two index loops. Both 

these graphs have two momentum loops. For our purposes both of these Feynman diagrams can also 
be generated by a simple pictorial algorithm: we splice together propagators of primitive diagrams 
using the vertex in Figure ||c, as displayed in Figure ^a and b. All graphs of the double-trace theory 
can be generated from primitive diagrams by this simple algorithm. Note that the number of index 
loops never changes when primitive diagrams are spliced by this pictorial algorithm. 

If a splicing of diagrams does not create a new momentum loop we say that the diagrams have 
been pasted together. This happens when the diagrams being spliced are originally disconnected 
as, for example, in Figure ||a. In fact because of momentum conservation, no momentum at all 
flows between pasted diagrams. If a new momentum loop is created we say that that the diagrams 
have been pinched. This happens when two propagators within an already connected diagram 
are spliced together as, for example, in Figure |2|b. In this example one momentum loop becomes 
two because momentum can flow through the double-trace vertex. Further examples of pinched 
diagrams are given in Figure |3| where the new loop arises from momentum flowing between the 
primitive diagrams via double-trace vertices. 
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(a) Paste 



(b) Pinch 



Figure 2: With the inclusion of the double-trace term we need new types of vertices. These can be 
obtained from the "primitive" diagrams associated with the pure single-trace superpotential by (a) pasting 
or (b) pinching. The vertices have been marked with a cross. 




Figure 3: More examples of "pinched" diagrams. 

To make the above statement more clear, let us provide some calculations. First, according to 
our operation, the number of double index loops never increases whether under pasting or pinching. 
Second, we can calculate the total number of independent momentum loops £ by £ = P — V + 1 
where P is the number of propagators and V, the number of vertices. If we connect two separate 
diagrams by pasting, we will have Ptot = (-Pi — 1) + (P2 — 1) + 4, Vtot = V% + V2 + 1 and 

£tot = Ptot - Vtot + l = h+£ 2 , (2.24) 

which means that the total number of momentum loops is just the sum of the individual ones. If 
we insert the double-trace vertex in a single connected diagram by pinching, we will have Ptot = 
P - 2 + 4, Vtot = V + I, and 

£tot = Ptot -V tot + l = £ + l, (2.25) 
which indicates the creation of one new momentum loop. 
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Having understood the structure of double-trace diagrams in this way, we can adapt the tech- 
niques of f|] to our case. The steps 1-6 as described in Sec. 2.1 go through without modification 
since they are independent of the details of the tree-level superpotential. However the steps 7-11 
are modified in various ways. First of all naive counting of powers of fermionic momenta as in step 
7 leads to the selection rule 

h > £, (2.26) 

where h is the total number of index loops and I is the total number of momentum loops. (The 
holomorphy and symmetry based arguments of || would lead to the same conclusion.) Since no 
momentum flows between pasted primitive diagrams it is clear that this selection rule would permit 
some of the primitive components to be non-planar. Likewise, both planar and some non-planar 
pinching diagrams are admitted. An example of a planar pinching diagram that can contribute 
according to this rule is Figure |2|b. However, we will show in the next subsection that more careful 
consideration of the structure of perturbative diagrams shows that only diagrams built by pasting 
planar primitive graphs give non-zero contributions to the effective superpotential. 

2.3 Which Diagrams Contribute: Selection Rules 

In order to explain which diagrams give non-zero contributions to the multi-trace superpotential 
it is useful to first give another perspective on the fermionic momentum integrations described in 
steps 7-9 above. A key step in the argument of @ was to split the glueball insertions up in terms of 
auxiliary fermionic variables associated with each of the momentum loops as in ( p.16 ). We will take 



a somewhat different approach. In the end we want to attach zero or two fields to each index 
loop, where p labels the index loop, and the total number of such fields must bring down enough 
fermionic momenta to soak up the corresponding integrations. On each oriented propagator, with 
momentum 7rj a , we have a left index line which we label pl and a right index line which we label 
Pr. Because of the commutator in (|2,14|), the contribution of this propagator will be 



exp(- Si (7T ia (^ PL) ~ W? m) ))- (2-27) 

Notice that we are omitting U(N) indices, which are simply replaced by the different index loop 
labels. In a standard planar diagram for a single-trace theory, we have one more index loop than 
momentum loop. So even in this case the choice of auxiliary variables in ( [2.27 ) is not quite the 



same as in ( 2.16| ), since the number of W a s is twice the number of index loops in ( |2.27| ) while the 



number of auxiliary variables is twice the number of momentum loops in (2.16). 

Now in order to soak up the fermionic ir integrations in ( |2.8| ), we must expand (2.27) in powers 
and extract terms of the form 

H,)K,)---K^ ( 2 - 28 ) 

where £ is the number of momentum loops and all the pi are distinct. The range of p is over 1, . . . , h, 
with h the number of index loops. In the integral over the anticommuting momenta, we have all h 
W(p) appearing. However, one linear combination, which is the 'center of mass' of the W( p ), does 
not appear. This can be seen from (2.27): if we add a constant to all W( p ) simultaneously, the 
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propagators do not change. Thus, without loss of generality, one can set the W( p ) corresponding 
to the outer loop in a planar diagram equal to zero. Let us assume this variable is WrM and later 
reinstate it. All W( p ) corresponding to inner index loops remain, leaving as many of these as there 
are momentum loops in a planar diagram. It is then straightforward to demonstrate that the W 



appearing in ( 2.16| ) in linear combinations reproduce the relations between propagator momenta 



and loop momenta. In other words, in this "gauge" where the W corresponding to the outer loop is 
zero, we recover the decomposition of W a in terms of auxiliary fermions associated to momentum 
loops that was used in j| and reviewed in ( [2.16 ) above. 



We can now reproduce the overall factors arising from the fermionic integrations in the planar 



diagrams contributing to ( 2.17 ). The result from the tt integrations is some constant times 



P =i 

Reinstating W^) by undoing the gauge choice, namely by shifting 

W(„) - W(p) + (2.30) 



for p = 1, . . . , h — 1, ( 2.29 ) becomes 

i 

n( w ( P ) +n w 2 - ( 2 - 31 ) 

P =i 

The terms on which each index loop there has either zero or two W insertions are easily extracted: 

k=l \p^k ) 

In this final result we should replace each of the W? \ by S, and therefore the final result is of the 
form 

hS h -\ (2.33) 



as derived in Q and reproduced in ( 2.17 ). 

Having reproduced the result for single-trace theories we can easily show that all non-planar 
and pinched contributions to the multi-trace effective superpotential vanish. Consider any diagram 
with t momentum loops and h index loops. By the same arguments as above, we attach some W( p ) 
to each index loop as in ( j2~27l) , and again, the 'center of mass' decouples due to the commutator 
nature of the propagator. Therefore, in the momentum integrals, only h — 1 inequivalent W( p ) 
appear. By doing £ momentum integrals, we generate a polynomial of order 2£ in the h — 1 
inequivalent W?U. This polynomial can by Fermi statistics only be non-zero if £ < h — 1: VV|L is 
zero for all p. Therefore, we reach the important conclusion that the total number of index loops 
must be larger than the number of momentum loops 

h > £ (2.34) 
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while the naive selection rule ( 2.26| ) says that it could be larger or equal. 



Consider pasting and pinching k primitive diagrams together, each with hi index loops and ti 
momentum loops. According to the rules set out in the previous subsection, the total number of 
index loops and the total number of momentum loops are given by: 

h = Y J K ; £>J2 £ i ( 2 - 35 ) 

i i 

with equality only when all the primitive diagrams are pasted together without additional momen- 
tum loops. Now the total number of independent Ws that appear in full diagram is Ylii^i ~ -0 
since in each primitive diagram the "center of mass" W will not appear. So the full diagram is 
non- vanishing only when 

£<^Oi-l). ( 2 - 36 ) 

i 

This inequality is already saturated by the momenta appearing in the primitive diagrams if they 
are planar. So we can conclude two things. First, only planar primitive diagrams appear in the 
full diagram. Second, only pasted diagrams are non-vanishing, since pinching introduces additional 
momentum loops which would violate this inequality. 

Summary: The only diagrams that contribute to the effective multi-trace superpotential are 
pastings of planar primitive diagrams. These are tree-like diagrams which string together double- 
trace vertices with "propagators" and "external legs" which are themselves primitive diagrams of 
the single-trace theory. Below we will explicitly evaluate such diagrams and raise the question of 
whether there is a generating functional for them. 

2.4 Summing Pasted Diagrams 

In the previous section we generalized steps 7 and 8 of the the single trace case in Sec. 2.1 to the 
double-trace theory, and found that the surviving diagrams consist of planar connected primitive 
vacuum graphs pasted together with double-trace vertices. Because of momentum conservation, no 
momentum can flow through the double-trace vertices in such graphs. Consequently the fermionic 
integrations and the proof of localization can be carried out separately for each primitive grapg, 
and the entire diagram evaluates to a product of the the primitive components times a suitable 
power of #2 j the double-trace coupling. 

Let Gi, i = 1, . . . , k be the planar primitive graphs that have been pasted together, each with 
hi index loops and l{ = hi — 1 momentum loops to make a double-trace diagram G. Then, using 



the result (2.18|) for the single-trace case, the Schwinger parameters in the bosonic and fermionic 



momentum integrations cancel giving a factor 

hoson ■ If ermion = H(NhiSf) = N k S^ h ^ J[ h U (2.37) 

i i 

where the last factor arises from the number of ways in which the glueballs S can be inserted into 
the propagators of each primitive diagram. Defining C(G) = Y\ i hi as the glueball symmetry factor, 
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k{G) as the number of primitive components, h{G) = ^ hi as the total number of index loops and 
i(G) = Yli h = h(G) — k(G) as the total number of momentum loops, we get 



hoson ■ Ifermion = H(NhiSj) = N h ^~ 1 ^ C '(G) . (2.38) 
i 

We can assemble this with the Veneziano-Yankielowicz contribution contribution for pure gauge 
theory [25] to write the complete glueball effective action as 

W eff = -NS(log(S/A 2 ) - 1) + C(G)^{G)N h(G ^- e ^S e(G \ (2.39) 

G 

where J~{G) is the combinatorial factor for generating the graph G from the Feynman diagrams 
of the double-trace theory. Notice that in our discussion, we have set 52 = m = 1> so A 2 in this 
equation is in fact mA 2 which matches the dimension of S. We can define a free energy related to 
above diagrams as 

To = j2 :F ( G ) sh(G) - ( 2 - 4 °) 

G 

J-q is a generating function for the diagrams that contribute to the effective superpotential, but 
does not include the combinatorial factors arising from the glueball insertions. In the single-trace 
case that combinatorial factor was simply Nh(G) and so we could write W e ff = N(dJ r o/dS). Here 
C{G) = \\hi is a product rather than a sum h(G) = ^2 hi, and so the effective superpotential 
cannot be written as a derivative of the free energy. 

Notice that if we rescale 52 to g2/N, there will be a N~ ( fe ( c? ) _1 ) factor from k(G) — 1 insertions 
of the double-trace vertex. This factor will change the N h( - G ^~ 1 ^ dependence in ( |2.3S| ) to just N 
for every diagram. This implies that the matrix diagrams contributing to the superpotential are 
exactly those that survive the large M limit of a bosonic U (M) Matrix model with a potential 

V($) = <7 2 Tr($ 2 ) + 54 Tr($ 4 ) + | 2 r Tr($ 2 )Tr($ 2 ). (2.41) 

In Sec. 4 we will compute the large M limit of a such a Matrix model and compute the free energy 
J-q in this way. 

Below we will compute this effective action ( |2.39D to the first few orders. In Sec. 3 we will show 
that it is reproduced by an analysis based on the Seiberg-Witten solution of N = 2 gauge theories. 
In the single-trace case Dijkgraaf and Vafa argued that the large- iV limit of an associated Matrix 
model carries out the sum in ( 2.39p , or equivalently, that the matrix model free energy provides a 



generating function for the perturbative series of matrix diagrams contributing to the exact field 
theory superpotential. In Sec. 4 we will show that the well known double-trace Matrix models that 
have large-iV limits do sum up the same "planar pasted diagrams" that we described above and 
give the free energy defined by ( 2.4C| ). However, unlike the single-trace case, the Matrix model will 



not reproduce the the combinatorial factors C{G) appearing in (p. 3 
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2.5 Pert ur bat ive Calculation 



Thus equipped, let us begin our explicit perturbation calculations. We shall tabulate all combi- 
natoric data of the pasting diagrams up to third order. Here C(G) = Y\ { hi and T(G) is obtained 
by counting the contractions of <J>s. For pure single-trace diagrams the values of J-(G) have been 



computed in Table 1 in [26 1, so we can utilize their results. 
2.5.1 First Order 

To first order in coupling constants, all primitive (diagram (b)) and pasting diagrams (diagram (a)) 
are presented in Figure |||. Let us illustrate by showing the computations for (a) . There is a total of 
four index loops and hence h = 4 for this diagram. Moreover, since it is composed of the pasting of 
two primitive diagrams each of which has h = 2; thus, we have C(G) =2x2 = 4. Finally, T = 52 
because there is only one contraction possible, viz, Tr( , ii , )Tr( , ii , ). 




(a) g : 




(b) 



Figure 4: All two-loop primitive and pasting diagrams. The vertices have been marked with a cross. 
In summary we have: 



diagram 


(a) 


(b) 


h 


4 


3 


C{G) 


4 


3 


HG) 


92 


2g 4 



(2.42) 



2.5.2 Second Order 

To second order in the coupling all primitive ((c) and (d)) and pasting diagrams ((a) and (b)) 
are drawn in Figure || and the combinatorics are summarized in table (2.43). Again, let us do 
an illustrative example. Take diagram (b), there are five index loops, so h = 5; more precisely it 
is composed of pasting a left primitive diagram with h = 3 and a right primitive with h = 2, so 



C{G) = 2x3 = 6. Now for F{G), we need contractions of the form Tr($$ $ $)Tr($ $)Tr ("¥¥"); 
there are 4 x 2 x 2 = 16 ways of doing so. Furthermore, for this even overall power in the coupling, 
we have a minus sign when expanding out the exponent. Therefore J~{G) = —16g294 for this 
diagram. 

In summary, we have: 



diagram 


(a) 


(b) 


(c) 


(d) 


h 


6 


5 


4 


4 


C{G) 


8 


6 


4 


4 


HG) 


-453 


-I69294 


-2gj 


-16gj 



(2.43) 
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Figure 5: All three-loop primitive and pasting diagrams. The vertices have been marked with a cross. 



2.5.3 Third Order 



Finally, the third order diagrams are drawn in Figure ||. The combinatorics are tabulated in ( 2.44 ). 
Here the demonstrative example is diagram (b), which is composed of pasting four diagrams, each 
with h = 2, thus h(G) = 4 x 2 = 8 and C(G) = 2 4 = 16. For F(G), first we have a factor ^ from 



i 1 „ i 1 1 1 

the exponential. Next we have contractions of the form Tr(<I><3?)' ;i Tr(<I> <I>)Tr(<I> $)Tr($ $); there are 

2 3 x 4 x 2 ways of doing this. Thus altogether we have T(G) = for this diagram. 
In summary: 



diagram 


(a) 


(b) 


(c) 


(d) 


(e) 


(/) 


(9) 


h 


8 


8 


7 


7 


7 


6 


6 


C(G) 


16 


16 


12 


12 


12 


8 


8 


HG) 




32^3 
3 92 


64^94 


329194 


649194 


1289291 


128^291 


diagram 


(h) 


W 


(J) 


(k) 


(0 


(m) 




h 


6 


6 


5 


5 


5 


5 




C{G) 


8 


9 


5 


5 


5 


5 




HG) 




649291 


1289! 


32 3 
3 9 A 


649| 


256 _3 

~9a 





(2.44) 



2.5.4 Obtaining the Effective Action 

Now to the highlight of our calculation. From tables ( 2.42j ), (2.43), and ( 2.44} ) we can readily 
compute the effective glueball superpotential and free energy. We do so by summing the factors, 
with the appropriate powers for S, in accordance with ( |2.39 , 2~40| ). 
We obtain, up to four-loop order, 

T(G)S h( -^ 



^0 



E 

G=all diagrams 



(2 9 4 + 92S)S 3 



2(99 4 2 + 89492S + 2g 2 2 S 2 )S 4 



16 



+ — (54gl + mglg 2 S + 30g 4 gtS 2 + 5g 3 2 S 3 )S 5 + • 



(2.45) 



and subsequently, 

W eff = -iV5(log(5/A 2 ) 



1)+ C{G)T(G)N h <®- l M&W 

G=all diagrams 
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Figure 6: All four-loop primitive and pasting diagrams. The vertices have been marked with a cross. 

= -NS(log(S/A 2 ) - 1) + (6 54 + 4g 2 N)NS 2 - (72gj + 96^452^ + 32g%N 2 )NS 3 
20 

+y(6 5 4 + 45 2 iV) 3 iV < S 4 + ---. (2.46) 

We will later see how this result may be reproduced from independent considerations, i.e., the 
effective action from the factorization of Seiberg-Witten curve and free energy from the correspond- 
ing matrix model. 
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2.6 Multiple Traces, Pasted Diagrams and Nonlocality 

Above we found that the only diagrams that contribute to the effective superpotential have zero 
momentum flowing though the double-trace vertex. Now observe that the double-trace term in the 
tree- level action can be written in momentum space as: 

V = J d 4 xTi{$ 2 (x))Tr{$ 2 (x)) 

= J d 4 Pl d 4 p2d 4 P3 d 4 p 4 Tr($(p 1 )<S>(p 2 ))Tr(<S>(p 3 )$( P ±))5(p 1 +p2+P3+p±). (2.47) 

Since no momentum flows through the double-trace vertices contributing to the superpotential, the 
delta function momentum constraint factorizes in our pasted diagrams as 

5 (pi +P2+P3+P4) ~ $(Pl +P2)5(P3 +Pa)- (2.48) 

Therefore for the purposes of computing the superpotential we might as well replace the double- 
trace term in the action by 



V = J d 4 Pl d 4 P2 d 4 P3 d 4 P4 Tr^(p 1 M P2 ))Trmp 3 M P4 ))5(p 1 + p 2 ) 5(p 3 +p A ) 

= j d 4 x d 4 yTr($ 2 (x))Tr(<S> 2 (y)). (2.49) 

The Feynman diagrams of this nonlocal theory include the ones that compute the superpotential 
in the double-trace, local theory. 

This fact suggests that correlation functions of chiral operators are position independent, as 
described in [||. Along with cluster decomposition, this position independence leads to the state- 
ment that correlators of operators in the chiral ring of an M = 1 theory factorize, which we use 
in the next section to write (Tr(<I> 2 ) 2 ) = (Tr(<I> 2 )) 2 . However it is subtle to establish the precise 
equivalence between factorization of chiral operators and the vanishing that we demonstrated of 
all except the pasted diagrams, which leads in turn to the nonlocal action in ( 2.49j ). We leave this 
potential connection for exploration in future work. 

3. The Field Theory Analysis 

In this section, we will show that in the confining vacuum the effective superpotential of the field 
theory discussed in the previous section is 

W eff = NA 2 + (6Ng 4 + AN 2 g 2 )A 4 . (3.1) 

After integrating in the glueball superfield and expanding the superpotential in a power series in 
S, (|3.1|) can be written 

20 

W ef f = -NS{log(S/A 2 ) - 1) + (6# 4 + 4g 2 N)NS 2 -2{6g 4 + Ag 2 N) 2 NS 3 + — (<og A + Ag 2 N) 3 NS 4 + • • • . 



(3.2) 
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We shall compare this expression for the low-energy gauge dynamics to the perturbative field theory 
computations in Sec. g. The two results, of course, are in concert. 

We begin by considering an J\f = 1 U(N) gauge theory with a single adjoint superfield <E> 
deformed from N = 2 by the tree-level superpotential (n < N) 



n+l 



Wtree = E 9rU r + ^g 2 ul, (3.3) 



r=l 

where 



u k := ±Tr($ k ). (3.4) 



The tree-level superpotential in (3.3) is more general than the one used in ( p. 23 ). Here, we allow 



single-trace terms at arbitrary powers of We shall specialize to the previous example at the end 
of our discussion. 

3.1 The Classical Vacua 

To find the classical vacua, we have to solve the D-term and F-term conditions. The D-term is 
proportional to Tr[<3?, $] 2 which is zero if is diagonal. Let the diagonal entries be Xi, i = 1, . . . , N. 
We still need to solve the F-term condition. In terms of the Xj, the tree-level superpotential is 

2 

(3.5) 



n+l N / N \ ■ 

r=l i=l \i=l / 



From this, the F-flatness condition reads: 

71+1 / N 

dx k 



= jr = E 9rX T k - 1 + ±92X k ( E X H ' k = 1 > N - (3 ' 6) 



r=l \i=l 



This is certainly different from the case without the double-trace term, where the F-term equations 
for different x k s decouple. Here, the eigenvalues interact with each other even at the classical level. 
To solve the (3.6), which may be recast as 



j n+l / N \ 

Xk r=l \i=l J 



(3.7) 



x k 

we take the RHS of Q 

/ N \ 

" (3.8) 



as an unknown constant for all N F-terms. This gives 

n+l 

^grx^ 1 + Cx k = Vfe. (3.9) 

r=l 

Now the F-terms are decoupled. We can solve this system just as we solve for the vacua of a field 
theory with only single-trace interactions p8[ simply by taking g 2 i— > g 2 + C. 
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As the F-terms are order n polynomials in x, we should generically expect n solutions for each 
eigenvalue Xk- The eigenvalues are the roots of the polynomial 

n+l n 

= ^grxl' 1 = g n+ i - ai). (3.10) 

r=l i=l 

If Ni of the eigenvalues are located at ai, where Yli Ni = N, the unbroken gauge symmetry is 

n 

U(N)^YlU(Ni). (3.11) 

i=i 

As ajS are a function of C, we need to impose the additional consistency condition that 

n 

4g 2 ^N?a* = C. (3.12) 

3=1 

To simplify the discussion, we henceforth focus on the special case where all of the XjS have 
the same value. The SU{N) part of the gauge group is unbroken, and will confine in the infrared. 

3.2 The Exact Superpotential in the Confining Vacuum 

We now proceed to find the exact superpotential in this confining vacuum |[27| . Let us recall the 
general philosophy of the method (see e.g., p9t , whose notations we adopt, for a recent discussion). 
A generic point in the moduli space of the U(N) M = 2 theory will be lifted by the addition of 
the the general superpotential (|3.3| ). The points which are not lifted are precisely where at least 
N — n mutually local monopoles become massless. This can be seen from the following argument. 
The gauge group in the M = 1 theory is broken down to n?=i U(Ni), and the SU(Ni) factors each 
confine. We expect condensation of Ni — 1 magnetic monopoles in each of these SU(Ni) factors 
and a total of N — n condensed magnetic monopoles. These monopoles condense at the points 
on the M = 2 moduli space where N — n mutually local monopoles become massless. These are 
precisely the points which are not lifted by addition of the superpotential. These considerations are 
equivalent to the requirement that the corresponding Seiberg-Witten curve has the factorization 

P N (x, uf - 4A 27V = H N ^ n (xfF 2n (x), (3.13) 

where Pn{x, u) is an order polynomial in x with coefficients determined by the (vevs of) the u^, 
A is an ultraviolet cut-off, and H and F are, respectively, order N — n and 2n polynomials in x. 

The N — n double roots place N — n conditions on the original variables Uk- We can parametrize 
all the (uk) by n independent variables ctj. In other words, the ays then correspond to massless 
fields in the low-energy effective theory. If we know the exact effective action for these fields, to 
find the vacua, we simply minimize S e ff. Furthermore, substituting (uk) back into the effective 
action gives the action for the vacua. 

Holomorphy and regularity of the superpotential as the couplings go to zero requires that there 
are no perturbative corrections to the tree-level superpotential. In addition, we assume that all 
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non-perturbative effects are captured in the Seiberg-Witten curve analysis discussed above. Then, 
we need to minimize 

n+l 

W exact = J2g r (u r ) + Ag 2 (ul). (3.14) 

r=l 



In general the factorization problem is hard to solve [30], but for the confining vacuum where all 
N — 1 monopoles have condensed, there is a general solution given by Chebyshev polynomials. 2 In 
our case, we have the solution 

AT W 

«p = - E CpC*^*-* (3.15) 

P q=0 

Notice that in ( p. 15 ), there is one free parameter z which is the field left upon condensation. Now 
we put it into the superpotential 

n+l 

W = J29rU r {z,A)+Ag 2 u 2 {z,A) 2 , (3.17) 

r=l 

solve and back-substitute z from dW/dz = to obtain the effective superpotential W e ff. Notice 
that in the above result, we have used 

{u 2 2 ) = (u 2 ) 2 . (3.18) 

This is true because u 2 is a chiral field, and cluster decomposition in the field theory lets us factor 
the correlation functions of operators in the chiral ring ||. 

Although the above procedure finds W e ff, it is not the best form to compare with our previous 
results because there is no gluino condensate S. To make the comparison, we need to "integrate 



in" [31] the glueball superfield as in |p9[ . 

The integrating in procedure is as follows (here we use the single-trace superpotential as an 
illustrative example of the technique). 
• We set A := A 2 , and use the equation 

aw n+1 [p/2] 

NS = A flA = E -C 2 v q Cl q z^^ (3.19) 

' r =2 q=l P 



to solve for A in terms of S. 
• Next, we find z by solving 



o w n+l [p/2] 

= IT = E^ E V —^-CfCl^-^bfi. (3.20) 



2 This was worked out first by Douglas and Shenker j32|, but here we use the results and nomenclature of Ferrari 

il. 
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Now the effective action for the glueball superfield S can be written as 



W s (S,g,A) = -Slog[ — ) +W tree (S,g,A), (3.21) 



which will reproduce the result 



^W S (S, A 2 ,g) = - hx(A/A 2 ) N . (3.22) 



3.3 An Explicit Example 

Let us work out the double-trace example that we are interested in solving. The superpotential is 

W = g 2 u 2 + 4,g 4 u 4 + Ag 2 u 2 2 , (3.23) 

(later, we can set g 2 = m = 1). Using 

u 2 = ^[ Z 2 + 2A 2 ], (3.24) 

n 4 = j[^ 4 + 12AV + 6A 4 ], (3.25) 



from ( |3.15| ), we obtain 



W = Nz 2 [^ + g 4 z 2 + 12 54 A 2 + 4Ng 2 A 2 + Ng 2 z 2 ] (3.26) 
+NA 2 [g 2 + 6g 4 A 2 + ^Ng 2 A 2 ]. 

From this we have following equations by setting A 2 = A: 

S = A[g 2 + (1204 + ±Ng 2 )z 2 + (12<? 4 + 8Ng 2 )A], (3.27) 

= — = z[g 2 + z 2 {Ag A + 4Ng 2 ) + 2A(12 54 + 4iV£ 2 )]. (3.28) 

We solve z = 0: 

A = -92 + y/gj + 45(12 g4 + 8Ng 2 ) 

2(12g 4 + 8Ng 2 ) ' 1 ' ' 

The effective action with S integrated in is 3 

W eff = Slo g (- + NA[g 2 + A(6g 4 + ANg 2 )}, (3.30) 

which will be the one used in comparison to our previous results. After minimizing this action, we 
find 

W z=0 = NA 2 [g 2 + 6g 4 A 2 + 4Ng 2 A 2 } , (3.31) 



Notice that in both formula (3.29) and (3.30), g$ and 32 combine together as 174 + =£-gn. So if we shift <?4 to 
94 + ^j-g2, the single-trace result will reproduce the double-trace result, and the effective action for the double-trace 
can be naively calculated from the DV prescription by partial differentiation of the glueball field S. 
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which is the promised result of (3T). Setting g 2 = 1, expanding ([3.29 ) in powers of S, and 
substituting into (|3.30| ), we get the second formula ( |3.2[) from the beginning of this section: 

20 

W ef f = -NS{log{S/A 2 ) - 1) + (654 + Ag 2 N)NS 2 - 2{Qg A + Ag 2 N) 2 NS 3 + — {6g 4 + Ag 2 N) 3 NS i + ■■■. 

(3.32) 

Crucial in matching the result of this calculation with the perturbative analysis in the previous 
section is the assumption of factorization (u 2 ) = {u2) 2 - This is equivalent to the vanishing the 
pinching diagrams in the perturbative analysis. 



4. The Matrix Model 

In Sec. 2 we demonstrated that the explicit field theory computation of the effective superpotential 
localizes to a certain sum of matrix diagrams. All of these diagrams are constructed by pasting 
planar single-trace diagrams together with double-trace vertices in such a way that no additional 
momentum loops are created. By examining the scaling of these diagrams with N in ( 2.3S| ) we also 



observed that these are precisely the diagrams that would survive the large-M limit of a U(M) 
bosonic Matrix model with a potential 

= g 2 Tr($ 2 ) + 54 Tr($ 4 ) + g:Tr($ 2 )Tr($ 2 ). (4.1) 

The extra factor of 1 /M multiplying g 2 , in comparison with the field theory tree- level superpotential 
(|1.1| ) is necessary for a well-defined 't Hooft large-M limit. This is because each trace, being a sum 
of eigenvalues, will give a term proportional to M. So to prevent the double-trace term from 
completely dominating the large-M limit we must divide by an extra factor of M. Fortunately, 



this a well known model and was solved more than a decade ago [22, 33]. Below we will review 
this solution, compare its results to our field theory calculations, and then generalize to other 
multi-trace deformations. 



4.1 The Mean-Field Method 

The basic observation, following [J22|j, that allows us to solve the double-trace matrix model ( [4,1| ). 
is that in the large-M limit the effects on a given matrix eigenvalue of all the other eigenvalues can 
be treated in a mean field approximation. Accordingly we compute the matrix model free energy 
T as 

exp(-M 2 ^) = J d M2 (cD) exp{-M(lTr(<D 2 ) + g 4 Tr(cD 4 ) + g 2 )}, (4 . 2) 

= / 1] dX, exp{M(-i £ A? - § £ Xf - A £ m + £ log ^ _ A . | } 

i i i i ijtj 

Here A are the M eigenvalues of $ and T is the free energy, which can be evaluated by saddle 
point approximation at the planar limit. The log term comes from the standard Vandermonde 
determinant. This matrix model is Hermitian with rank M in the notation of |?5| (of course 
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as mentioned earlier, we should really consider GL(M, C) matrices though the techniques hold 



equally). We have introduced an extra factor of M in the exponent on the right hand side of (4.2) 



by rescaling the fields and couplings in ( [4.1[ ) in accordance with the conventions of [22]. 
The density of eigenvalues 

1 M 

: = M E <H A - A ») ( 4 - 3 ) 

i=i 

becomes continuous in an interval (—2a, 2a) when M goes to infinity in the planar limit for some 
a S M + . Here the interval is symmetric around zero since our model is an even function. The 
normalization condition for eigenvalue density is 

r-2a 

d\p(X) = 1 . (4.4) 

-2a 

We can rewrite (4.2) in terms of the eigenvalue density in the continuum limit as 

/ M f-2a i 

Y[d\iexp{-M 2 ( / dAp(A)(-A 2 + o 4 A 4 ) (4.5) 
i=l J -2a ^ 

p2a r2a p2a 

+ g 2 ( dXp(X)X 2 ) 2 - / / dXdpp(X)p(p)ln\X-p\)} . 



J -2a J — 2a J — 2a 

Then the saddle point equation is 

1\ + 204 A 3 + 2g 2 cX = P ^ dp-^, (4.6) 

2 J -2a A - /i 

where c is the second moment 

/2a 
dXp(X)X 2 (4.7) 
-2a 

and P means principal value integration. 

The effect of the double-trace is to modify the coefficient of A in the saddle point equation. We 
can determine the number c self-consistently by (4.7). The solution of p(X) to (4.6) can be obtained 



by standard matrix model techniques by introducing a resolvent. The answer is 



p{X) = i Q + 2g 2 c + Ag 4 a 2 + 2g 4 X 2 ^j Aa 2 - X. (4i 



Plugging the solution into (4.4) and (4.7) we obtain two equations that determine the parameters 
a and c: 

a 2 (l + 4g 2 c) = 1 - 12o 4 a 4 , (4.9) 
160452a 8 + (12<74 + 4ff 2 )a 4 + a 2 - 1 = 0. (4.10) 
Substituting these expressions into ([4.5|) gives us the free energy in the planar limit M — > 00 



as: 

r2a 1 r f2a 



/la 1 p pZa 

dXp(X)(-X 2 + C/4A 4 ) +g 2 c 2 - / / dXdpp(X)p(p)ln\X- p\. 
-2a 2 J J- 2a 

(4.11) 
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One obtains 



^(<74, h) - HO, 0) = \{a 2 - 1) + (6 54 a 4 + a 2 - 2)<? 4 a 4 - ~ log(a 2 ) . (4.12) 



Equation ( |4. 12 ) together with ( 4.1 0| ) give the planar free energy. We can also expand the free 



energy in powers of the couplings, by using ( |4.10| ) to solve for a 2 perturbatively 

a 2 = 1 - (12,74 + 452) + (28%f + 176 5 452 + 32g ) (4.13) 
-(8640 5 | + 7488gjg 2 + 2496g 4 g 2 2 + 320^1) + • • • . 



Plugging this back into (4.12) we find the free energy as a perturbative series 



1 6 

Fo = T{g 4 , 52) - ^(0, 0) = 2g 4 + g 2 - 2{9gj + 8g 4 g 2 + 2g 2 ) + — (54 5 | + 66gjg 2 + 30 94 5l + 5$) + ■ ■ ■ ■ 

(4.14) 

Comparing with ( 2.45| ) we see that T reproduces the explicit computation of the generating 



function of "planar pasted" field theory diagrams in Sec. 2.5. In matching the two we have to restore 
the proper powers of the glueball S into ( (4.14|) . First recall that to keep the relevant diagrams in 
the matrix model, we have inserted i to the double-trace term in (4.1). Therefore to compare 



with ( 2,45| ), we need to rescale g 2 in ( [4.14 ) to g 2 M = g 2 S where we have effectively identified the 



glueball S in the field theory with M in the matrix model. In addition, we should re-insert powers 
of M into ( 4.14 ) by loop counting. The first two terms in ( [4.14j ) have three index loops so we need 



to multiply them by M 3 = S 3 . The third term has four index loops and fourth term, five and 
hence we respectively need factors of S 4 and S 5 . With these factors correctly placed into ( 4.14| ), 
we recover ( [2.45 ) completely. 



This verifies our claim that the diagrams surviving the large-M limit of the matrix model fl4.1|) 
are precisely the graphs that contribute to effective action of the field theory with the tree-level 
superpotential (1.1). Nevertheless, as we have already discussed in Sec. 2.4 we cannot compute the 



effective superpotential of the field theory W e ff(S) by taking a derivative dJ^o/dS, because the 
combinatorial factors will not agree. In this way the double-trace theories differ in a significant 
way from the single-trace models discussed by Dijkgraaf and Vafa |J. We can pose the challenge 
of finding the field theories who effective superpotential is computed by the Matrix model ( |4.12j ). 

4.2 Generalized Multi- Trace Deformations 

In fact the mean field techniques of the previous subsection can be generalized to solve the general 
multi-trace model. Below we illustrate this by solving the general quartic Matrix model; as discussed 
above, it is an interesting challenge to find a find a field theory whose effective superpotential these 
models compute. 

Specifically, let us consider the Lagrangian 

C = <7 2 Tr($) 2 + /u(Tr($)) 2 + ^(Tr$) 2 (Tr($ 2 )) + ^ 2 (Tr$) 4 + 2^ 3 (Tr$)(Tr($ 3 )), (4.15) 

which exhausts all quartic interactions. 
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The one-matrix model partition function 



Z M = J [DX] exp j-M 2 (l - J dx dy log |A(x) - A(y)|^ | , (4.16) 
gives the saddle point equation 

f 26 u(t) 

4g 4 \ 3 + 6v 3 c 1 \ 2 + 2(g2 + 2g2C2 + v 1 c 2 1 )\ + 2(fic 1 + v l c 1 C2 + 2v 2 cl + V3C 3 ) = 2P / dr-r-^-, (4.17) 

J -2a X~T 



where the moments c& are defined as 

/•2b 

c k = dTu(r)T k , 

J -2a 

c = / otu(t) = 1. (4.18) 

J -2a 

Note that we have introduced the separate upper and lower cut parameters a and b as opposed to 
the standard symmetric treatment because u(X) is not of explicit parity (such asymmetric examples 
have also been considered in ^6|). When a = b one can recast fl4.17|) into a Fredholm integral 
equation of the first kind and Cauchy type, which affords a general solution as follows [pCl] 



1 f a u(t)dt . . , , 1 f a a 2 - 1 2 \ v(t)dt C ri . 

p - ■7 J - = -v(x)=>u(x) = - / \ —r -jy I t^— + -=== (4.19) 



7rJ_ a t-x n J_ a \\ a 2 - x 2 } t-x Va 2 - x 2 

for some constant C. When a ^ b we can use the ansatz: 

1 



u(A) = - (^A 2 + 5A + C) v/ (2a + X)(2b - A) (4.20) 

7T 



with the constants matching the coefficients in the LHS of ( 4.17 ) as 

2A = 2g 4 , (4.21) 
2a A- 2Ab + 2B = 3u 3 ci, 
-a 2 A-2aAb- Ab 2 + 2aB -2bB + 2C = g 2 + 2g 2 c 2 + v\c\, 
a 3 A + a 2 Ab — a Ab 2 — Ab 3 — a 2 B — 2ab B — b 2 B + 2aC — 2bC = \ic\ + vic x c 2 + 2v 2 c\ + u 3 c 3 . 

We see a well-behaved u(X) which is zero at the end-points and vanishes outside the support 
(-2a, 26). 

We now need to check the consistency of our mean-field method. This simply means the follow- 
ing. Considering the definition of Cj in ( [4. 18 ), the definitions ( 4.21| ) actually constitute a system of 
equations for A, B, C, a, b because each Cj on the RHS, through ( |4.18| ), depend on A, B, C, a, b. To 
( 4.21] ) we must append one more normalization condition, that Co = Jj^a u {X)dX = 1. Therefore we 
have five equations in five variables which will fix our parameters in terms of the seven couplings. 
Our mean-field method is therefore self-consistent. It would be interesting to find a role for such 
exactly solvable models in the physics of four-dimensional field theories. 
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5. Linearizing Traces: How To Identify the Glueball? 



In previous sections we showed that the field theory computation of the effective superpotential of 
a double-trace theory as a function of the glueball S localized to summing Matrix diagrams. In 
the end, the only the "pasted" diagrams that contributed, namely certain tree-like graphs obtained 
by pasting together planar single-trace graphs with double-trace vertices in such a way that no 
momentum flows through the latter. After verifying the result via the N = 2 Seiberg-Witten 



solution, we demonstrated that the sum of these diagrams given as a series in ( 2.39 ) is not computed 
by the large M limit of a U(M) Matrix as one would have naturally hoped. Nevertheless, we may 
wonder if there is some Matrix model that sums the series of pasted diagrams. In this section we 
take up the challenge of finding such a Matrix model. 

Since the authors 0] and || have proven that the superpotential of a single-trace gauge theory 
can be computed from an associated Matrix model, we seek to construct our double-trace theory 
from another single-trace model. Recall that we are considering the tree-level superpotential 

W tree = ^ 2 Tr(cD 2 ) + g 4 Tr(<S> 4 ) + g 2 (Tr(d> 2 )) 2 . (5.1) 

Now consider another theory with an additional gauge singlet field A 

W tree = i( 52 + 45 2 A)Tr($ 2 ) + 54 Tr($ 4 ) - g 2 A 2 . (5.2) 

It is easy to see that integrating out A in ( |5.2j ), which amounts to solving d ^g% ee = and back- 
substituting, produces the double-trace theory ( |5.1| ). 

The advantage of (|5.2[) is that it consists purely of single-trace operators. The first two terms 
will generate an effective potential Wgingie^j 5), as a function of A and the glueball superfield S 
(the subscript "single" refers to the fact that this is the superpotential for the model without the 
double-trace term and with an A dependent mass term). Then 

W cS (A, S) = W singlc (A, S) - g 2 A 2 . (5.3) 

The exact superpotential for the glueball superfield S for the double-trace theory then follows by 
integrating A out, i.e. solving dw ^ la _ 2g 2 A = for A and substituting in (|5.3|) . Since single- 
trace theories are directly related to Matrix model we might hope to use this construction with an 
added auxiliary field A to find an auxiliary Matrix model that sums the pasted diagrams of the 
double-trace theory. 

5.1 Field Theory Computation of W sing i e (A, S) and Pasted Matrix Diagrams 

In this section 4 we will discuss how the superpotential for the double-trace theory can be computed 
in field theory from the linearized model (|5.2j ). First, observe that the superpotential for an adjoint 
theory with an additional gauge singlet ( |5.2| ) localizes to summing matrix integrals, since the 
arguments of Q that are reviewed in Sec. |2| go through essentially unchanged. To compute the 

4 We thank Cumrun Vafa and Ken Intriligator for communications concerning the material in this section. 
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effective potential as a function of both A and the glueball S, we need to sum superspace Feynman 
diagrams with insertions of both A and W a , with both of these treated as background fields. Since 
we are only interested in contributions to the superpotential, we can restrict ourselves to constant 
background A. Then it is easy to verify that the entire analysis in Sec. ^ goes through for the 
theory fl5.2|) , with the double-trace coupling g 2 set to zero and a shift in the mass of the field 
<]?, viz g 2 — ► 92 + 4^2- In particular, the computation of the effective superpotential W s i ng \ e (A, S) 
localizes to summing matrix diagrams and there is some free energy Single in terms of which 
W single = NdT singlc (S,A)/dS. 

Let us verify that this procedure will yield the correct double-trace result when we integrate A 
out. Making the g 2 — > g 2 + 4^2 with 52 = in the known single-trace result Q3.30D and ( |3.29| ), we 
find the effective superpotential 

A 2 

W singlc (A, S) = NSlog(-) + NA(( 92 + Ag 2 A) + 6g 4 A), (5.4) 

where A is determined by the quadratic equation 

12 54 A 2 + ( 52 + 4g 2 A)A = S. (5.5) 

Now we can integrate A out and obtain the superpotential for the glueball superfield S. We solve 
dWes/dA = dW sing i c /dA - 2g 2 A = for A. This is a simple calculation from (0), (|J): 

= (92 + 4<?2^4 + 12 54 A - + 4g 2 NA 

= ig 2 2NA, 

where in last step we have used (|5.5|) . So the solution to dW s i ng \ e /dA — 2g 2 A = is 

A = 2NA. (5.7) 



Plugging A = 2NA into (p.4|) , (|5.5|) and (|5.3|) , we find the effective glueball superpotential to be 



W eS = NS dog(^) + iVA( 5 2 + 452^A + 654A), (5.8) 
with A determined by the quadratic equation 

(1204 + 8g 2 N)A 2 + g 2 A = S. (5.9) 
Equations (5.8) and (|5.9[) are of course the double-trace effective glueball superpotential we com- 



puted previously in ( 3.29j ) and ( |3.30|) . 



Why does this procedure reproduce precisely the sum of pasted diagrams that contribute to 
the double-trace superpotential in ( 2.39| )? From the point of view of perturbation theory, we are 



doing the following. If we treat A as a constant, we should simply sum the planar diagrams in the 
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theory with a quartic superpotential, and after doing so, we obtain the superpotential J d 4 xd 2 9 W 
with 

W = ^connected planar 

(S,g 2 + 4g 2 A, gi ) -92 A 2 . (5.10) 

Next, we should integrate out A. To do so, we write A = Aq + A, where Aq solves dW/dA = 0. 
We see that W becomes 

W = ^connected planar 

(S, g 2 + 4g 2 A ,g 4 ) - g 2 A 2 + c 2 A 2 + c s A 3 + .... (5.11) 

What is the meaning of integrating over A? From the diagrammatic point of view, A is the field 
that allows momentum to flow through the g 2 vertices. All diagrams where such momentum flow is 
prohibited are taken into account by the background value Aq. Thus, picking Aq takes the pasting 
process into account, whereas the further integrals over A should correspond to pinching diagrams. 
We already know that these latter diagrams should vanish from our diagrammatic analysis, and 
therefore we should simply drop all terms involving A. The final answer for W is thus 

W = W conri e C ted planar 

(S,g 2 + 4g 2 Ao,g 4 ) - g 2 A 2 . (5.12) 

One can also see directly that integrating out A gives no contribution to the superpotential. In the 
diagrams that one can write down, there will be many loops of A, but there are no vertices that 
can absorb any fermionic momentum, and therefore these diagrams do not yield any contribution 
to the superpotential. 

It is an interesting exercise to verify explicitly that ( |5.12f ) is a generating diagram for pasted 
diagrams (which are all tree graphs), made out of building blocks corresponding to W con nected planar ■ 

5.2 Matrix model perspective 

Above we argued that the methods of show that above linearization of the double-trace defor- 
mation via introduction of an auxiliary singlet field A leads to a theory whose superpotential be 
computed by a matrix model. 

First observe that the double-trace matrix model partition function can be linearized in traces 
by the introduction of an auxiliary parameter A, over which we integrate: 

Z = exp(-M 2 ^oubi e) = J ^ Wexp{ _ M( l g2Tr( ^ 2) +g4Tr( ^ 4) + ~ (Try)) 2 )} ({ug) 

= JdA d M \$) exp{-M(i( 52 + 4g 2 A)Ti{<& 2 ) + 94 Tr($ 4 ) - Mg 2 A 2 )}. 

This is is Matrix model analog of the statement that the double-trace field theory can be generated 
by integrating out a gauge singlet. In terms of the free energy of the single-trace matrix model, 
this can be written as 

exp(-M 2 ^ doublc ) = J dAexp(-M 2 J single + M 2 g 2 A 2 ). (5.14) 

Hence to obtain the free energy of the double-trace matrix model, we need to solve the equation 
ggq^Ag) _ 2g 2 A = for A and substitute in jrdoubic = ^ingie^^ _ where we have 
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used the identification from |J that S ~ M. The resulting expression for the double-trace matrix 
model is, of course, the same as that obtained by mean field methods in section (|4.l| ). However, as 
emphasized in the previous sections, the derivative of this free energy with respect to S does not 
yield the correct superpotential for the field theory. 

Let us now contrast this with a different matrix model construction suggested by the field 
theory analysis in the previous subsection. Consider the Matrix partition function 

Z = exp(-M 2 ^ singlc ) = J d M \$)e W {-M(±(g 2 + 4g 2 A)Tr(<!> 2 )+g 4 Tr(<S> 4 )}, (5.15) 

where A is now treated as a fixed parameter of the Matrix model in analogy with the constant A 
appearing in the field theory superpotential. As we explained above, the arguments of 0] applied 
to the linearized model fl5.2|) show that in terms of J^ msle (A, S) (with S ~ M as in ||) 



a ^-single 

W eS (A, S) = -NS(log(S/A 2 ) - 1) + N—l- 



dS 



-g 2 A 2 , (5.16) 

constantA 



where W e g is the field theory superpotential in (|5.3| ). Note that we have not integrated out A at this 
stage. Since the single-trace matrix model has reproduced W e ff (-^i S), it is manifest that integrating 
out A will correctly produce the double-trace superpotential as a function of the glueball, just as 
it did in the field theory. 

Identifying the glueball: Something remarkable appears to have happened here. If we start 
with the single-trace Matrix model (5.13) and integrate out the singlet A, we find the free energy of 



the double-trace model as indicated in Q5.14 ), and dJ-Q° nble /dS does not reproduce the field theory 



superpotential. However, if we first differentiate with respect to S and then integrate out A we 
reproduce the field theory result. Of course, the direct field theory computation described above 
indicates to us that this is the right order in which to do things. But notice that the difference 
between the two procedures essentially amounts to correctly identifying the glueball. Apparently to 
identify S in the Matrix model we must linearize the theory and then use the prescription given in 
H]. Of course, differentiating the single-trace free energy with respect to S at constant A translates 
in the double-trace theory into some complicated operation which would in effect identify the field 
theory glueball in that Matrix model. But it is challenging to identify what this operation is. 

5.3 General Multi-trace Operators 

Finally, we show that the procedure of introducing auxiliary parameters to linearize traces in a 
double-trace theory can be extended to a general multi-trace model. Consider the term 

Tr($ mi )Tr($ m2 ). (5.17) 

in the superpotential. We can rewrite this in terms of single trace terms by introducing four gauge 
singlet fields Ai , i = 1 • ■ • 4 as 

W 2 = 3 (a\ +A\ + A 1 A 2 + AiTr$ mi + A 2 Tr(<J? m2 + -^=A 3 TT^ mi -A\ + A 4 Tr$ m2 - A^j . (5.18) 
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Integrating out Ai by setting dWz/dAi = 0, solving for Ai and substituting in ( 5.18| ) yields the 
double trace superpotential fl5.17[ ). 
To generate a term of the form 

Tr($ mi )Tr($ m2 )Tr($ m3 ) . (5.19) 

we iterate the above procedure twice, i.e. we introduce additional gauge singlet fields Bi, i = 1 • • • 4, 
and consider the theory with a superpotential 



W 3 = 3(B( + Bi + B t B 2 + B^r^ + B 2 W 2 + — J B 3 Tr$ m3 - B$ + ^=B 4 W 2 -Bi). (5.20) 



where W2 is defined in ( |5.18 ). Integrating out A{ and Bi for i = 1---4 yields the term ( 5.19j ). 



Generalization to terms with more traces is obvious. 



6. Conclusion 

We have studied an N = 1 U(N) gauge theory with adjoint chiral matter and a double-trace 
tree-level superpotential. We found by direct computation that the computation of the effective 
superpotential as a function of the glueball superfield localizes to summing a set of matrix integrals. 
The associated set of Matrix diagrams have the structure of tree diagrams in which double-trace 
vertices are strung together by "propagators" and "external" legs are that themselves connected 
single-trace diagrams. We showed that the Seiberg-Witten solution to J\f = 2 field theories computes 
an effective superpotential for the double-trace theory that matches our direct analysis. The use 
of factorization in our Seiberg-Witten analysis, namely that (Tr(<I> 2 ) 2 ) = (Tr(<I> 2 )) 2 , was equivalent 
in our perturbative computations to the vanishing of any diagrams where extra momentum loops 
were introduced by the double-trace vertices. Next, we showed that the large-M limit of the 
standard double-trace U(M) Matrix model does sum up the same set of matrix diagrams, but the 
combinatorial factors are different from those appearing in the field theory. In particular, the field 
theory superpotential is not computed by a derivative of the matrix model free energy as in ||. Put 
another way, a simple manipulation of the free energy of the standard double-trace matrix model 
does not give a generating function for the field theory superpotential. Finally we demonstrated 
how a multi-trace field theory can be linearized in traces by the introduction of auxiliary singlet 
fields. We showed the associated Matrix model, which is linear in traces also, computes the field 
theory superpotential as a function of both the glueball and the new singlets. The basic subtlety, 
then, lies in the correct identification of the field theory glueball as a variable in a Matrix model. 
Our results raise several challenges: 

1. We found that while a multi-trace Matrix model did not directly compute the superpotential 
of a multi-trace field theory, we could sum the necessary diagrams by introducing auxiliary 
singlet fields in an associated single-trace model. The basic subtlety involves correct identifi- 
cation of the field theory glueball in the Matrix model. How is this done in general, and what 
is the underlying principle driving the identification? 
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2. 



Does the large- N limit of the standard double-trace Matrix model compute the superpotential 
for some M = 1 field theory? 



3. We expect that all our results generalize easily to multi-trace theories — it would be nice to 
check this. 

4. We have worked in the vacuum with an unbroken gauge group. It would be good to generalize 
our arguments to the other vacua with partially broken gauge symmetry. 

5. In Sec. 2.6 we pointed out that there is an intriguing connection between the contributions 
made by multi-trace vertices to the superpotential of a local theory and certain Feynman 
diagrams of an associated nonlocal theory. It would be very interesting to flesh out this 
connection. 

6. In an U(N) theory with adjoint the operator Tt{^ K ) with K > N decomposes into a sum 
of multi-trace operators. This decomposition can receive quantum corrections as discussed in 
H|. How do our arguments generalize to this case? 

7. We can also add baryon-like operators like det(3>) to the superpotential (for theories with 
fundamental matter in the context of matrix models, baryons were studied by ||, 0). Such 



operators also decompose into sums of products of traces, and are very interesting because, 
even without fundamental matter, they can give rise to an open string sector in Yang-Mills 
theory as opposed to the standard closed string sector that the 't Hooft expansion leads 



us to expect [37, Bq]. It would be useful to understand in this case how and whether the 



computation of holomorphic data in such a theory localizes to sums of Matrix integrals. 

In addition to these directions there are some interesting applications that arise from known 
facts about the large-A^ of the standard double-trace U(N) Matrix model. This theory is related 
to two-dimensional gravity with a positive cosmological constant and displays phase transitions 
between branched polymer and smooth phases of two-dimensional gravity p2j] . Presumably such 
phase transitions manifest themselves as interesting phenomena in a four-dimensional field theory. 
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